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Abstract 

We report on the results of a computer search for primes p which divide 
an Harmonic number with small N > 1. 
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1 Introduction 

Before its complete proof by Andrew Wiles, a major result for the first case 
of Fermat’s last theorem (FLT), that is, the assertion of the impossibility of 
_l_ yP = 2 :^ in integers x,y,z, none of which is divisible by a prime p > 2, 
was the 1909 theorem of Wieferich [25] that the exponent p must satisfy the 
congruence 


2P-1 _ 1 

qp{2) := - = 0 (mod p), 

P 

where qp{2) is known as the Fermat quotient of p, base 2. This celebrated result 
was to be generalized in many directions. One of these was the extension of 
the congruence to bases other than 2, the first such step being the proof of an 
analogous theorem for the base 3 by Mirimanoff m in 1910. Another grew 
from the recognition that some of these criteria could be framed in terms of 
certain special Harmonic numbers of the form 

[P/N] ^ 

^ip/Ni := X 7 

for TV > 1, and with [-J denoting the greatest-integer function. There are nice 
historical overviews of these developments in and [H] (pp. 155-59), which 
go into greater detail than we attempt here. 
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1.1 Congruences for l^[p/Arj 

At the time when Wieferich’s and Mirimanoff’s results appeared, it was already 
known that three of these Harmonic numbers had close connections to the Fer¬ 
mat quotient, and satisfied the following congruences (all modulo p): 


H[p/2] — 2 • qp{2) ( 2 ) 

^Lp/3J ~ ~ 2 

l^[p/4j = • qp{‘2). (4) 

All these results are due to Glaisher; those for N = 2 and 4 will be found in 
m, pp. 21-22, 23), and that for A^ = 3 in ([ID], p. 50). Although it was not 
the next Harmonic number criterion published, it will be convenient to dispense 
next with the case N = 6. The apparatus needed to evaluate this Harmonic 
number appears in a 1905 paper of Lerch m, p. 476, equations 14 and 15), 
but the implications of Lerch’s result were long overlooked, and only realized 
in 1938 by Emma Lehmer m, pp. 356fF), who gave the following congruence 
mod p: 


^[p/6] — 2 • 9p(2) 2 • 5^(3). (5) 

The fact that the vanishing of H^p/^j mod p is a necessary condition for the 
failure of the first case of FLT for the exponent p is an immediate consequence 
of the theorems of Wieferich and Mirimanoff. In the present study, we use 
Lehmer’s congruence in the equivalent form 

H[p/6i =-2' ( 6 ) 

obtained from ([^ by applying in reverse the logarithmetic and factorization 
rules for the Fermat quotient given by Eisenstein |7]. This expression reveals 
that divisors p of L7[p/6j are instances of the vanishing of the Fermat quotient 
mod p for composite bases, a problem which has notably been studied in the 
ongoing work of Richard Fischer [^. 

The four congruences above exhaust the cases of i7[p/Arj that can be evalu¬ 
ated solely in terms of Fermat quotients. Historically, the case N = 5 also has 
its origins in this era, and was in fact introduced before that of TV = 6. In 1914, 
Vandiver m proved that the vanishing of both i7[p/5j and qp{5) are necessary 
conditions for p to be an exception to the first case of FLT. Unlike the four cases 
already considered, here the connection of the Harmonic number with FLT was 
discovered before any evaluation of it (beyond the definitional one) was known. 
It was almost eighty years later that the connection between these results would 
become apparent. The ingredients needed for the evaluation of this Harmonic 
number were presented in a 1991 paper by Williams ([DD], p. 440), and almost 
simultaneously by Z.H. Sun m, pt. 3, Theorems 3.1 and 3.2); and though 
they do not write out the formula explicitly, it is clearly implied to be 
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( 7 ) 


5 5 

77Lp/5J = -4 • 9p(5) - 4 • Pp-il,)/P (mod p), 

where is the Fibonacci quotient (OEIS A092330), with F a Fibonacci 

number and 3- Jacobi symbol. In light of Vandiver’s theorems on F[^p/^j and 

qp{5), this result immediately established that the vanishing of the Fibonacci 
quotient mod p was yet another criterion for the failure of the first case of FLT 
for the exponent p. The fact was announced almost immediately in the cele¬ 
brated paper by the Sun brothers m, which gave fresh impetus to an already 
vast literature on the Fibonacci quotient; and in its honor the primes p which 
divide their Fibonacci quotient were named Wall-Sun-Sun-Primes. These re¬ 
main hypothetical, as not a single instance has been found despite tests to high 
limits [TI?] . 

The remaining known formulae for the type of special Harmonic numbers 
in which we are interested are of much more recent origin. In some cases they 
were discovered simultaneously, or nearly so, by more than one researcher; and 
we hope we have not done injustice to any of the participants. Apart from a 
few published formulae which are apparently in error or under determined, we 
have the following congruences (all modulo p) which are undoubtably correct: 


HLp/ 8J =-4-gp(2)-2.C/p_(|)(2,-l)/p (8) 

Hyp/w\ = -2 • qp{2) - ^ • qpib) - ^ • Fp_^i.^/p (9) 

i7Lp/i2j = -3 • qp{2) - ^ • qp{3) - 3 • (I) • t7p_(|)(4, l)/p (10) 

JJ[p/24J = — 4 • qp{2) — - ■ qp{3) — 4 • C/p_(2)(2, —l)/p 

- 3 ■ (I) ■ i)/? - 6 - (p) q-(|)(io. i)/p 

where the are Jacobi symbols, and 

• is as before the Fibonacci quotient (OEIS A092330) 

• ?7p_(2^(2, —l)/p is the Pell quotient (OEIS A000129) 

• Vp_(3)(4, l)/p is a quotient derived from the Lucas sequence 1, 4, 15, 56, 
209, ^ . (OEIS A001353) 

• ?7p_(6)(10, l)/p is a quotient derived from the Lucas sequence 1, 10, 99, 
980, 9701, ... (OEIS A004189) 
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• (I) = (-IjO'’-'"* 

• (f) = 

, ^6^ ^ (_l)L(P+5)/12j^ 

The result for N = 8 is derived from the 1991 paper by Williams ([53]) P- 
440), with an equivalent result also appearing in Sun ([55], pt. 3, Theorem 3.3). 
To the best of our knowledge, at the time of the discovery of this formula in 
1991 the vanishing of modulo p was not recognized as a condition for the 

failure of the first case of FLT for the exponent p, and this only became evident 
with the appearance of the 1995 paper of Dilcher and Skula [3] discussed below. 

The result for TV = 10 is due to a 1992 paper by Z.H. Sun (|1S|, pt. 3, 
Theorem 3.1). Since the vanishing of all the individual components was by then 
known to be a necessary condition for the failure of the first case of FLT for the 
exponent p, the same thing was immediately seen to be true for H^p/iQj. 

The result for = 12 is also derived from the 1991 paper of Williams ( [29]. 
p. 440), and that for TV = 24 is from a 2011 paper of Kuzumaki & Urbanowicz 
([Il],p. 139). 

In their landmark joint paper of 1995, Dilcher and Skula |3| proved among 
other things that the vanishing modulo p of was a necessary criterion 

for the failure of the first case of FLT for the exponent p, for all N from 2 to 46. 
This result gives retrospective interest to the evaluations of and Lf[p/i 2 j i 

and furnishes the main motivation for the present study. 

1.2 The case of Af = 9 

As discussed in section]^ below, when calculating the sort of special Harmonic 
number with which we are concerned, it is important to avoid by any means 
possible the calculation of the long runs of modular inverses suggested by the 
definition (Q. In the early days of this study Tl[p/ 9 j was obtained from H^p/i^j , 
requiring the calculation of inverses for ^ ~ ^ of th® numbers in the range 

{l,p}. This is more efficient than obtaining il[p/i8j from Lf[p/ 24 j, the nearest 
neighbor that can be calculated by formula, which requires the calculation of 
inverses for ~ ^ ^ of the numbers in the range {l,p}. However, assum¬ 

ing 71 |^p/i8j will be calculated anyway, lf[p/ 9 j can be derived from it without 
calculation of any additional inverses, and with negligible cost. 

The device proposed is somewhat convoluted, but worthwhile. In 1991, Zhi- 
Hong Sun (|5^, Theorem 2.7) gave an evaluation of 

^ k ’ 

k=l 

which is in fact trivially equivalent to 
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Lp/9J 

»( 1 . 18 )= E 

fc=Lp/18j+l 

to use a well-known notation for this sum. Consider the recurrence defined by 


/(0)=l,/(l) = 0,/(2) = 2, 

f{n) =^f{n-2) - f{n-3) (n > 3). 

Sun showed that there is a composition, which we shall designate Z{p), defined 
by 

Z(p) = s(l, 18) - 2gp(2) (modp), (13) 

such that 


Zip) = 


/(p+l)-2 

P 

fip) - /(P - 1) - 2 


-fiP + 1) - fip) + /(P - 1) - 2 


if p = ±1 (mod 9) 
if p = ±2 (mod 9) 
if p = ±4 (mod 9). 


(14) 


The evaluation of the sequence fin) can be performed efficiently in PARI using 
the matrix exponentiation ([0,3, —1; 1, —1,1; 1,0, l]^n)[l, 1], with the calcula¬ 
tions being performed modulo p^ throughout. Thus, noting that by definition. 


we have from (131 that 


Hip/9\ = 7?[p/i8j + s(l, 18), 


Hip/g] = -ff[p/i8j + Zip) -k 2gp(2) (mod p). (15) 

No direct evaluations of iP[p/i8j and ^^[p/gj comparable to that for s(l, 18) are 
known, but Sun’s result allows efficient evaluation modulo p of either of these 
special Harmonic numbers with respect to the other. Compared with a scheme 
in which the two quantities are calculated separately, the derivation of iPj^p/gj 
from requires the calculation of | as many modular inverses, while 

the derivation of iP|^p/i8j from if[p/ 9 j requires the calculation of only | as many 
inverses. However, we choose the former path because we also require the values 
of i7|^p/7vj with 18 < fV < 24. 
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2 The Divisibility of Harmonic Numbers 

It may be helpful to distinguish the purpose of the present study with that of the 
more general problem of the divisibility of Harmonic numbers. Eswarathasan 
and Levine [5] note that all primes greater than 3 divide the Harmonic numbers 
of indices p — 1, p{p — 1), and — 1, and Wolstenholme’s theorem states that 
for a prime p > 3, p^|iLp_i. These results may be inverted to give three general 
rules for the divisibility of Harmonic numbers: 

1. Hn is divisible by (n + 1)^ if n + 1 is a prime > 3; 

2. Hn is divisible by if the latter is a prime > 3; 

3. Hn is divisible by if the latter is a prime > 3. 

Eswarathasan and Levine define harmonic primes as primes that divide only 
the three aforementioned Harmonic numbers (OEIS A092101), and anharmonic 
primes as those that divide additional Harmonic numbers (OEIS A092I02). 
From the fact that we consider only cases Q where p|iL|^p/jvj with N > 1, 
it will be evident that any such p is anharmonic, and that we are seeking a 
subset of anharmonic primes that divide Harmonic numbers of relatively small 
index; for example, from Table we see that N = 2 is solved by p = 1093, 
implying that 1093|iL546, while V = 46 is solved by p = 11731, implying that 
11731|i7255- In contrast, the work of Boyd [T] and Rogers m entails, in part, 
finding Harmonic numbers of large index divisible by relatively small primes, 
such as the case Il|i7i0ii84977i8552i49i29684042i7247- The overlap between our 
results and theirs is slight. 

Indeed, our results entail only a minority of the anharmonic primes, because 
they require that p divide an Harmonic number i7„ with n < p — 1, and thus 
belong to a subset of the anharmonic primes which has sometimes been called 
the Harmonic irregular primes (see OEIS A092194 and the Wikipedia entry for 
“Regular prime”). Although the standard definition of such a prime p is that 
it divide an Harmonic number with n < p — 1, it may be noted that any 
prime satisfying this criterion must in fact divide some Harmonic number 
of smaller index n < (p — l)/2, since by symmetry Hp_i_n = mod p. The 
inequality sign in this condition is strict because if n = (p — l)/2, then p is 
a Wieferich prime and by ([^ and (4) must likewise divide the smaller 
But this inequality also gives a sharp upper bound on the index of the least 
Harmonic number divisible by an Harmonic irregular prime p, since among the 
first 999 primes we have n= (p — 3)/2 for p = 29, 37, 3373 (see OEIS AI25854). 
In general, there does not seem to be any way of recognizing an Harmonic 
irregular prime other than by exhaustively testing it as a divisor of Harmonic 
numbers with n < (p — 3)/2 (though one would of course not attempt to do 
this by actual division except for very small p). 

In turn, our results entail only a minority of the Harmonic irregular primes, 
because we must have m = [p/N\ for some N. Nonetheless, for some small 
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primes this relationship is satisfied by more than one value of N; for example, 
with p = 137, p\H^, and [_p/23j = [p/24j = [p/25j = [p/26j = [p/27j = 5. 

Taking into account the stringency of the conditions on p, it is perhaps 
unsurprising that few instances have been found where a given p divides distinct 
Harmonic numbers of the special type under consideration. By and (4) all 
Wieferich primes p divide both i7[p/4j and H^p/ 2 \, but otherwise we have found 
only two instances where a non-Wieferich prime p < 72,000,000 divides two 
distinct Harmonic numbers of index m in such a way that m = lp/N\ is satisfied 
for some N < 1000 in each case: for p = 761, p divides both iJg (with N = 85 
through 95) and H 23 (with N = 32, 33), and for p = 845921, p divides both 
77ioii (with N = 836) and i7i524 (with N = 555). Perhaps further solutions 
exist with larger N, but we expect such primes to be rare. 

There is, however, a sense in which our results do relate to the general prob¬ 
lem of the divisibility of Harmonic numbers. While the formulae (|^ through 
(11) are framed in terms of p, they may conversely be seen as divisibility con¬ 
ditions on Hn- For example, with N = 2 (the Wieferich primes), we seek cases 
where the numerator of is divisible by a prime 2n-|-l, 4n-|-l, or 4n-|-3. The 
still unresolved case N = 5 asks whether it is possible for Hn to be divisible 
by a prime of one of the forms 5n + 1, 5n + 2, 5n -|- 3, or 5n -|- 4, and the still 
unresolved case N = 12 whether it is possible for to be divisible by a prime 
of one of the forms 12n -|- 1, 12n -|- 5, 12 ?t, -|- 7, or 12n -1-11. Incidentally, we also 
tested the numerators of Hn for divisibility by any of these linear forms without 
the restriction that the divisors be prime, and found only one additional solu¬ 
tion for n < 10, 000, for a divisor of the form 12n -|- 1: Hiq is divisible by 121, 
which is of course 11^. For the same ranges of p as covered by Table below, 
the following table shows, for small k, all linear forms for which there exists no 
known p = kn + r {r < k) dividing an Harmonic number Hn- 
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Table 1: Linear forms for which no known prime p = kn + r divides an 
Harmonic number (n < 24); starred rows contain all possible 
forms of p 


k 

r 

* 5 

1, 

2, 

3, 

4 

7 

1, 

3, 

4, 

5, 

8 

3, 

7 



9 

1, 

5, 

7, 

8 

10 

1, 

3 



11 

1, 

2, 

3, 

4, 

* 12 

1, 

5, 

7, 

11 

13 

1, 

2, 

3, 

4, 

14 

1, 

3, 

5, 

9 

15 

1, 

2, 

4, 

8, 

16 

5, 

7, 

9, 

11 

* 17 

1, 

2, 

3, 

4, 

* 18 

1, 

5, 

7, 

11 

19 

1, 

2, 

3, 

4, 

* 20 

1, 

3, 

7, 

9, 

21 

1, 

4, 

5, 

8, 

22 

1, 

3, 

7, 

9, 

23 

1, 

2, 

3, 

4, 

24 

11 

, 13, 

23 


The example noted above of 121 dividing Hiq underscores the fact that if fc — 1 
is a prime, then (fc — l)^|iLfc _2 and {k — 1)^ = 1 mod k. Thus, although there 
are no known prime divisors of of such forms as 12n + 1, 14n + 1, 18n + 1, 
and 20n + 1, it is not because such divisors are algebraically impossible. 


3 Computational Considerations 

In the cases where such formulae as (|^ through exist for the kinds of spe¬ 
cial Harmonic numbers with which we are concerned, it would be difficult to 
overemphasize their superiority for computational purposes over na’ive applica¬ 
tion of the definition of the Harmonic number as a sum of reciprocals Q . As we 
noted in [I] , the limits attainable by the two methods differ by several orders of 
magnitude. A striking example is furnished by the cases of (calculated 

by a particular simple formula) versus (a smaller sum by definition, but 

necessarily calculated by brute force). Of the two, the calculations for i?[p/6j 
were begun more than a year later, since (as will be explained below) they 
overlapped previous work which there initially seemed no purpose in repeating. 
However, the calculations for i?[p/6j now extend more than 50,000 times as far 
as those for i?[p/ 7 j. 




The nine cases of ^^[p/Arj represented by formulae have all been tested to 
high limits (Table without finding solutions for ^f[p/ 5 j and and 

this suggests the question of whether the latter two cases ever vanish modulo p. 
These cases have resisted tests to limits that are at least four orders of magnitude 
greater than the least zeros (wherever one is known) of with N < 46. 

These limits are at least two orders of magnitude greater than the least zeros 
(whenever one is known) of all the Fermat quotients qp{b) with base b < 100, 
than all three known zeros of the Pell quotient (see OEIS A238736 for its zeros), 
and than the lesser of the known zeros of the Lucas quotient C/p_(3)(4, l)/p (see 
OEIS A238490 for its zeros). But if there are mathematical reasons for the 
observed non-vanishing of 7L[p/5j and iL[p/i 2 j, they remain elusive. 

On the other hand, we are not inclined to attach as much significance to the 
fact that there are still no known zeros of for N = 17,18,20,29,31,43, 

as these cases require brute-force processing and have still not been particularly 
well tested. Explicit computation of Harmonic numbers, even using modu¬ 
lar arithmetic, is an unappealing scenario; as noted by Schwindt, for fixed N 
it has algorithmetic complexity of order p log p. Where such computation is 
unavoidable, it can be somewhat accelerated by employing the elementary re¬ 
arrangement 




3 

2 


Ln/2J 

E 


i=L«/4J+i 


1 

j 



fe=Ln/2j-|-l 
k odd 


which reduces the number of terms by about |. The term Hin/A\ could likewise 
be decomposed in a similar manner, though for the sake of simplicity this was not 
done in the present study; and more sophisticated rearrangements are possible, 
where all the even terms are successively generated from the odd terms. But 
such devices can bring about only modest improvements in runtime. 

As to the solution of the remaining cases with N = 17,18, 20, 29, 31,43, there 
are indeed heuristic reasons for doubting the efficacy of the present means of 
testing. If one considers the sequence of least divisors p of Harmonic numbers 
H^p/pjj for N = 2 through 46 (see Table |^, arranged in ascending order of p and 
displayed on a logarithmetic scale (see Figure 1), it is readily seen that the points 
tend to describe a line of positive curvature. Thus, within the limits studied, it 
has been found that the sequence of least divisors p grows at a superexponential 
rate. This portion of our data, accumulating since January 2014, now represents 
about 2 years of processing time on a typical desktop computer with a 3.20 GHz 
processor. Perhaps a few more such p can be found using the present methods 
and equipment, but it does not seem reasonable to expect that all six will be 
(assuming they exist). 
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Figure 1: Least divisors p < 319, 900,000 of Harmonic numbers iL^p/ivj for 
N = 2 through 46, arranged in ascending order of p; horizontal 
axis = order of discovery; vertical axis = log p. 


4 Previous Calculations 

The following account does not aim to be complete, but we hope that it includes 
all significant contributions to the problem. In referring to earlier work, for the 
sake of brevity we focus mainly on the results of successful searches. In |3], 
we noted that a statement in pp. 389-390 of m, implying that = 0 

mod p has a solution p < 2,000 for every N between 2 and 46 other than 5, is 
incorrect (at least if we require p > N to ensure that is not a vacuous 

sum). Indeed, as reported below, for N = 17,18,20,29,31,43 there are no 
solutions with p < 319, 900, 000, and for iV = 5,12 there are no solutions with 
p < 9,520,000,000,000. 

The divisors p of 77[p/2j and 77[p/4j are the Wieferich primes (OEIS A001220). 
The first of these, 1093, was found in 1913 by Meissner, and the second, 3511, in 
1922 by Beeger. The Wieferich primes have inspired some of the most intensive 
numerical searches ever conducted, but no further instances have been found. 
The current search record doubtless belongs to the ongoing test by PrimeGrid 
PO] , which had reached p < 484,357,700, 000, 000,000 as of 6 October 2015. 

The divisors p of are the Mirimanoff primes (OEIS A014127). The 

first of these, 11, was announced by Jacobi in 1828 m, the second, 1006003, by 
Kloss in 1965 [13]. No further instances have been found. The current search 
record appears to be p < 970,453, 984,500, 000, held by Dorais and Klyve |B]. 

We have found no computational literature on i7[p/5j apart from the negative 
result in [3| that it does not vanish modulo p for p < 2,000, though inevitably, 
reports of negative results are more difficult to locate. One wonders why such 
results were not reported by Schwindt [24] . since he could have obtained them 
as a byproduct of his study of the interval 77[p/5j — ^iny case, we 

have tested 77[p/5j to a much greater limit than that attainable by Schwindt’s 
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method, without finding any instance where it vanishes modulo p. 

The divisors of ^f[p/6j (OEIS A238201, but ignoring solutions with p < 6) 
were investigated in 1983 by Schwindt [23], whose study employed a rather na’ive 
method and only reached p < 600,000, finding the single nontrivial case p = 61. 
We know of no further computations relating to these numbers until they were 
obtained in another guise as the zeros of gp(432) by Richard Fischer [5], whose 
ongoing test, which had reached p < 65,442,059,485,993 as of 30 September 
2015, has found two further solutions. The vast range studied by Fischer has 
not been completely retested, and initially we carried our calculations just far 
enough to rediscover his second solution. The three solutions are included in 
Tables]^ andbelow. At a rather late stage in the present study we decided to 
continue the replication of his calculations, for the sole reason that (to the best 
of our knowledge) they have been performed only once. (A similar study by 
Lezak HZ], conducted in apparent unawareness of Fischer’s work, extends only 
top < 35,000,000,000.) 


5 The Present Calculations 

The present study, which supersedes the corresponding section of ]S|, does not 
attempt to extend existing calculations in the cases N = 2,3,4, 6, and relies 
on known values for those cases in the tables below. Although Cikanek ]2] also 
gives similar conditions for the first case of FLT for N up to 94 (with certain 
qualifications), it was decided to focus computing resources on the original range 
considered by Dilcher and Skula, with N ranging from 2 to 46. 

The search for solutions of p|i/[p/jvj has now been carried to 319,900,000 
for all values of N in the range, and solutions have been found for all but eight 
cases out of 45 (see Table below). For values obtainable by the formulae Q 
through (111, the work was performed in PARI, using matrix exponentiation 
with modular arithmetic to evaluate the required recurrence sequences, and 
it has now been run almost continuously on a typical desktop computer since 
January 2014, except in the case of = 6 which was only begun in May 2015. 
The following search limits have been attained: 


Table 2: Search limits on divisors p of Harmonic numbers H [p/jvj for which 
a formula exists 


N 

limit of p 

5 

9,520,500,000,000 

6 

14,015,000,000,000 

8 

2,335,500,000,000 

10 

2,335,500,000,000 

12 

9,669,000,000,000 

24 

2,335,500,000,000 
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For those values of N necessitating a brute-force search, the work was per¬ 
formed concurrently on a second machine, which has likewise been running 
from January 2014 to the present, with minor interruptions. After testing sev¬ 
eral programming languages, the best runtimes were again found to be achieved 
in PARI, which is very efficient at calculating modular inverses. To expedite 
the computations, values of were generally obtained either by addition 

or by subtraction from the nearest neighbor that could be calculated by for¬ 
mula. Thus, JJ[p/ 7 j was obtained from from iF[p/i 2 j, and all 

cases with > 12 from FJ|^p/ 24 j. Initially (as reported in an earlier version 
of this study) H^p/g^ was obtained from H^p/ig^, but later the more efficient 
strategy discussed in section [F^ above was implemented. The calculations were 
separated into multiple concurrent runs; hence the limits so far attained vary, 
as follows: 

Table 3: Search limits on divisors p of Harmonic numbers H [p/]vj for which 
no known formula exists 


A 

limit of p 

7 

333,000,000 

11 

349,000,000 

9; 13-23 

319,900,000 

25-46 

333,330,000 


In all, solutions remain unknown for A = 5,12,17,18,20,29,31,43, and so for 
as long as feasible the work will continue. 

The results of the searches follow. First (Table |^, we report separately 
those results for which the search-limits are given in Table above, since these 
have particular interest and in some cases coincide with OEIS sequences. Then 
(Table we report all results combined, including those in Table We believe 
that the zero p = 31251349243 for A = 24 is the largest presently known for 
any value of A > 1. 
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Table 4: Divisors p of Harmonic numbers for which a formula exists 


N 

P 

OEIS reference 

2 

1093, 3511 

A001220 

3 

11, 1006003 

A014127 

4 

C 

1093, 3511 

A001220 

0 

6 

61, 1680023, 7308036881 

A238201 

8 

269, 8573, 1300709, 11740973, 241078561 


10 

227, 17539, 4750159 


12 

— 


24 

137, 577, 247421, 307639, 366019, 
5262591617, 31251349243 



Table 5: Divisors p of Harmonic numbers for all N, 2 through 46 


N 

P 

N 

P 

2 

1093, 3511 

25 

137 

3 

11, 1006003 

26 

137, 67939 

4 

1093, 3511 

27 

137, 23669 

5 

— 

28 

20101 

6 

61, 1680023, 7308036881 

29 

— 

7 

652913 

30 

27089407 

8 

269, 8573, 1300709, 

31 

— 


11740973, 241078561 

32 

761 

9 

677, 6691 

33 

761 

10 

227, 17539, 4750159 

34 

1553 

11 

246277, 1156457 

35 

4139, 4481, 4598569 

12 

— 

36 

1297 

13 

43214711 

37 

1439, 26833 

14 

2267, 6898819 

38 

2473, 3527, 4047089 

15 

134227 

39 

407893 

16 

38723, 38993, 4292543 

40 

509, 177553 

17 

— 

41 

509, 151883 

18 

— 

42 

509, 190657 

19 

521, 911 

43 

— 

20 

— 

44 

6967, 27361 

21 

1423, 5693, 5782639, 212084723 

45 

609221 

22 

2843 

46 

11731 

23 

137, 264391 



24 

137, 577, 247421, 307639, 366019, 




5262591617, 31251349243 
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Table 6: Least divisors p < 319, 900,000 of Harmonic numbers -H[p/jvj for 
N = 2 through 46, arranged in ascending order of p 


p 

N 

11 

3 

61 

6 

137 

23 

137 

24 

137 

25 

137 

26 

137 

27 

227 

10 

269 

8 

509 

40 

509 

41 

509 

42 

521 

19 

677 

9 

761 

32 

761 

33 

1093 

2 

1093 

4 

1297 

36 

1423 

21 

1439 

37 

1553 

34 

2267 

14 

2473 

38 

2843 

22 

4139 

35 

6967 

44 

11731 

46 

20101 

28 

38723 

16 

134227 

15 

246277 

11 

407893 

39 

609221 

45 

652913 

7 

27089407 

30 

43214711 

13 
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